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A dielectric medium consisting of roughly polarized molecules is treated as a 3D 
disordered spin system (spin glass). A microscopic approach for the study of sta- 
tistical properties of this system on micrometer space scale and nanosecond time 
scale of standing electromagnetic wave is developed. Using ergodic hypothesis the 
initial 3D spin problem is reduced to two separate ID problems along external field 
propagation. The first problem describes the disordered spin chain system while the 
second one describes a disordered A^-particle quantum system with relaxation in the 
framework of Langevin-Schrodinger (L-Sch) type equation. Statistical properties of 
both systems are investigated in detail. Basing on these constructions, the coefficient 
of polar izability, related to collective orientational effects, is calculated. Clausius- 
Mossotti formula for dielectric constant is generalized. For dielectric permittivity 
function generalized equation is found taking into account Clausius-Mossotti gener- 
alized formula. 

PACS numbers: M 05.50.+q, 61.43.Fs, 61.46. +w, 71.23.Cq, 75.10.Hk, 75.10.Nr, 75.lG.Pq, 
77.22.Ej, 77.22. Ch 

I. INTRODUCTION 

The formation and control of periodic nanostructures in various type media (Media with 
Periodically Modulated Refractive Index (MPMRI)) are one of most important applied 
problems. First of all it is related to the possibility of the creation of compact UV or X- 



ray Free- 
example 



ectron Lasers (FEL) based on the principle of transition radiation (TR) (see for 
jj). Currently the following two problems are discussed intensively: 



1. A gas-plasma medium with periodically varied ionization density 



2. A specially periodical solid-state superlattice- 
posed of layers with different refracted indices 



ike (SSL) structures, which are com 
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TR is generated due to the difference in frequency-dependent dielectric constants (dielectric 
permittivity functions) of adjacent layers (recall that radiation power is proportional to 
[ei(u7) — €2(10)]'^). Therefore, an important problem is the possibility of the control of this 
difference by means of some external field. In other words, a problem is to construct a 
superlattice with dielectric constant difference between neighboring domains having a form 
[ei{w,g) — e2(w,g)]^, where g describes controlling parameters, ei{w,g) and e2(w,g) are 
dielectric permittivity functions in the neighboring regions. According to theoretical and 
experimental studies the periodical structures can be created in condensed matter by means 
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of external electromagnetic or acoustic fields (see in particular 

This idea have been applied recently for experimental TR generation j2^. In particular, 
it was shown that the electron beam of 20Mev, passing through amorphous silicon dioxide 
a—Si02 with the standing electromagnetic wave (of 10 GHz frequency and 3 fim wavelength) 
inside, produces anomaly high radiation. Preliminary investigations explain its high intensity 
as a result of the multiple passing of the electron beam through interfaces between regions 
with different dielectric constants. The theoretical study explained the appearance of one- 
dimensional superlattice order in random media by media polarization, caused by orientation 
relaxation of elastic molecular dipoles in external field |25l | . 

The main goal of this paper is systematic investigation of the relaxation processes and 
critical effects in the disordered 3D spin system type of a — Si02 under the influence of 
external electromagnetic field, which forms a standing wave in the medium. In particular, 
we develop a mathematical approach for description of statistical behavior of disordered 
3-D spin system along the direction of standing wave propagation on space-time scale of its 
period. 

In the framework of this approach we generalize Clausius-Mossotti equation for dielectric 
constant and the equation for frequency-dependent dielectric permittivity on the scale of 
space-time period of standing electromagnetic wave. 
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II. FORMULATION OF THE PROBLEM. NEW IDEAS AND BASIC 

FORMULAS 

Let us begin by discussing the fundamental problem of the dielectric constant space-time 
modulation (generation of dielectric constant superlattice with controlling parameters) in 
some types of amorphus materials. 

In this subsection, we shall give mathematical formulations for dielectric constant in 
the presence of an standing low electrical field. A particular attention will be pointed to 
investigation of dispersion properties of dielectric permittivity function where generation of 
collective orientation effects is possible under the standing low electromagnetic field. 

The starting point in our discussion will be Clausius-Mossotti relation for dielectric con- 
stant. It is known that in isotropic mediums (as well as in the crystals with cubic symmetry) 
the dielectric constant is well described by the Clausius-Mossotti formula (see for example 



^-^ = —yNlal, (1) 

m 

where A'^ is the concentration of particles (electrons, atoms, ions, molecules (or dipoles)) 
with the given m types of polarizabilities and correspondingly are coefficients of polar- 
izabilities. From this formula it follows that the static dielectric constant depends on the 
polarizabilities properties of the particles as well as on their topological order. In the exter- 
nal alternating field the homogeneity and the isotropy of the medium is often lost. Then, it 
is expected that the formula (0) will be applicable after the minor generalization. 

The object of our investigation will be solid state dielectrics t ype of amorphous silicon 
dioxide a — Si02- According to numerical ah initio simulations [2^, the structure of this 
type compound is well described by 3-D random network (see Fig. 1.) 

The white and black nodes in this figure correspond to different atoms while links between 
them correspond to covalent bonds. The redistribution of charges in outer electronic shells 
takes place because of the asymmetry of the bounded atoms. As a result some atoms 
acquire positive charge while others acquire negative ones. Thus this type of compound can 
be treated as a disordered 3-D rigid dipoles system (see Fig. 2). Below this system will be 
called a 3-D disordered spin system. For the description of amorphus media we use 3-D lattice 
with the lattice constant doiT) = {mo/po{T)}^^^, where mo is the molecule mass, po is the 
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FIG. 1: The structure of amorphous silicon dioxide (a — Si02) is described by 3D random network 
with covalent bonds. Every silicon vertex (gold sphere) has 4 edges and every oxygen vertex (red 
sphere) has 2 edges. 

density and T is the temperature. The lattice contains one random spin per elementary 
cell. Note that it has random direction as well as random location inside the cell. Suppose 




FIG. 2: The amorphous silicon dioxide (see Fig. 1.) maybe represented as a 3D lattice where in 
every cells randomly put one rigid dipole (spin). 

now that with the help of external electromagnetic filed a standing wave is formed in the 
medium: 



E{x] EQ,fl,ipQ) = 2EQsm{ipo)cos{kx), (po = Qto, 



(2) 
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where y^o and to are correspondingly the initial phase and time, Q is the wave frequency, 
k = 2'it/Xs and Xg is the wavelength (see Fig. 3). The following natural question arises. How 




FIG. 3: a) Disordered 3D medium under the standing electromagnetic wave. 

does dielectric constant change on scale of wavelength period and on time interval At <^ 
~ 10~^sec when the relaxation time of molecular dipoles is r ~ 10"^^ ^ 10~^^sec <^ 
This question is important because faster processes, like the transition or Cherenkov 
radiation, can go in the media. Note, that the time of relativistic electron is passing the 
wavelength (A^ ~ 10~^cm) of standing wave and the time of formation of transition or 
Cherenkov photons in this layer is smaller than 10~^^sec. This time interval is essentially 
lower than the time during which the standing wave is stationary. Since the wavelength is 
supposed to be much bigger than the inter-dipole distance ^ do, the Clausius-Mossotti 
relation still remains true. In this case the main problem is to calculate the polarizabilitiy 
coefficient related to orientational effects. 

Taking into account the external field, one can express the polarization of the matter at 
arbitrary point as the macroscopic self-consistent relation: 



P{r,do{T)) = ^p(r-f) = J2\j2nmam{l'-^EioM-f^], 1 = l{Uy,h), (3) 



where I is the vector of 3D lattice, p correspondingly the dipole moment of molecule. The 
second equation in Q contributes in the value of the dipole moment (spin). Note, that the 
number of the carriers of given polarization type in elementary cell is rim ~ {do{T))~^, am 
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coefficients of corresponding types polarizabilities taking into account the external field and 
Eioc is local field, the effective field, at the site of an individual molecule that causes the 
induced polarization. Each effect adds linearly to the net dipole moment per molecule which 
is a fact verified by experiments. Under the action of external field different polarization 
types arise in media. However as the simple analysis shows the value of polarizabilities 
coefficient determined by the orientational effects is essentially higher than others. 

Note that the elastic orientational polarizability coefficient in amorphous media adip{l — r) 
randomly depends on the cell location. This fact is the consequence of random orientation of 
the local field strengths Eiodl—r) with respect to the external field E{x; Eq, Q, ifo). Therefore 
all the terms in the right side of are basically known and well studied in literature (see 



for example j2f 



22l 



28| ) except from those which are connected with orientational effects. 
The orientational effects have a collective nature and are characterized by average value 
of random sum '^fadip{l — r) (sum of random coefficients of orientational polarizabilities). 
Multiplying both sides of the relation on the external field, we obtain: 

P(f, g)E(x,g) = -5t/(f,g) = 5^n^[5^a^(r-f)4oc(r-f)]E(x;Eo,fi,</Po), (4) 

m f 

where —5U{r, g) describes the potential energy of amorphous matter in the external field, 
symbol g shows parameters of standing wave [Eq, f2, (fo) (controlling parameters). Later we 
will consider the statistical properties of the medium in the direction of wave propagation. 

Taking into account formula (0}, one can obtain the following expression for the part of 
potential energy of the 3D spin system, which is related with orientation effects in external 
field: 

-SUdip{r, g) = ^dip{l - r)Eioc{l - r)E{x] Eq, n, v^o)- (5) 

r 

Let us pick out a layer with volume V = x Ly x Lz in the infinite crystal lattice, where 
Lx ~ (As) ^ do{T) and {Ly, L^) (oo, oo). It is easy to see that this volume is ffiled with 
infinite number of L2;-site random steric spin chains ( see Fig. 4). 

An important problem is now to calculation the mean value of the interaction potential 
between the spin layer and external field. Formally for it the following expression may be 
written: 

-5[/v(r,g) = -^5f/^^(/"l|f,g), 
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FIG. 4: The steric ID random spin chain system. 



-5ULAl±\r,g) = '^adipil - r)Eioc{l - r)E{x] Eo,n, v?o), l± = h{lx, ly), (6) 

la: 

where -SUL^{l±\f,g) is the interaction potential between the steric spin chain and external 
field. First, we take the mean value of the potential — 5f/y (r, g) on (y, z) plane: 

lim ± / 6Uyix,g)dS^ = hm ^ [ J^^U^Ahlx, g) dS^ = (5f/z.,(x; g))(i), 

b± ►oo ^_|_ J b± »oo ^_|_ J ^ — ' 

l± 

where S± = Ly x and (•••)(|) is averaging over all possible stable steric ID spin config- 
urations. The integral corresponds to the potential energy density -{6Ul^{x] g))(|) on {y, z) 
plane depends on the distance x. Taking into account the fact that the distribution of spins 
in the plan (y, z) is random but isotropic (see Fig. 5) it is simple to prove that in the limit 
of S'_L ^ oo it becomes the full self-averaging of spin system. 
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FIG. 5: The projection of 3D disordered spin system on the plan {y,z). 
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This implies that one may use Birgoff ergodic hypothesis [30| and averages by plan {y, z) 
interaction potential - f^L^; (2^; g) ) (|) obtains also by the integration of chain's intrinsic energy 
distribution: 

= Jlz(E;gME • 

where -6Ul^(Ei\x, g) shows the interaction potential energy between some steric spin chain 
with energy of E and external field and Z(E; g) denotes the energy distribution function 
(partition function) of ID steric spin chain configurations. The definition of distribution 
function will be given in section 4. Recall, that in ((7j) we take into account only negative 
values of E, because only for these values the spin chins maybe stable. 

Taking the average value of expression ((Zj), we obtain the following relation for the mean 
value of the potential increment: 

-{6Uv{r;g))v = -(5f/i^(E|x,g))(^.|) = adip{E^{x; g))^, 

(5f/i,(E|x,g))(,,i) = ^ r\6UL^{E\x,g))a^dx, (8) 

-l^x Jo 

where the bracket {■■.)x means the integration over x on a scale L^, the mean value of a full 
averaging potential energy of interaction between the spin chains and the external field is 
defined by term -((5f/i^(E|x, g))(^,j). 

Definition 1. In the relation (jHJ parameter ach will be named a collective polarization 
coefficient of steric ID spin chain. It is given by formula: 

_ (^^..(E|.,g))(.,) 

Note, that a^h is complex amount and describes the average value of steric spin chain 
polarizability with taking into account the lattice relaxation. When E'q it is simple to 
show that ach 0. 

Now one can write down the expression for the sum in the right part of the relation 
which takes into account the spin chains orientation effects in the external field: 

Nmam = N'nO^m + N^HOlch = Y " ^d.piV.' \g,, ^ ^ (10) 

mm m {E^ {x; Eq, Q, ^o)) ^ 
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where N^h = ndipN'^ = {do{T))~^N~^ -concentration of steric spin chains. From (fTUj) it 
follows that in areas where the field strength is small the orientation correction vanishes. 
In other words the external field on wavelength scale Q creates the alternating inhomo- 
geneities with different dielectric constants. These layers are stable in nanosecond scale 
At ~ 10-^°sec = 0.1ns. 

Using (fTUj) we can generalize Clausius-Mossotti equation on the space-time scale of stand- 
ing wave taking into account orientation effects: 



est{ 



A(g), where A(g) 



Air 



"a" 



(<5t/i^(E|x,g)(,,i) 



^st[g) + -^ - - -6 dlN, {E^{x;Eo,n,ipo)), 

Note that e^t is the label of stationary dielectric constant. 

Now we turn to the study of the dispersion property of frequency-depending dielectric 
constant (dielectric permittivity function). 

In the theory of dielectric relaxation, one writes the frequency-dependent dielectric con- 
stant e{uj) with Williams- Watts js^ function of dielectric relaxation by relation |32|: 
e(uj) - eoo , . , . 



g[a,uj) 



-iujt 



dF^{t)/dt 



;i2) 



f s Coo 

where eoo = e(i^ oo) is the high-frequency limit of the dielectric constant and = e{uj — *• 
0), the static dielectric constant which can be defined from generalized Clausius-Mossotti 
equation. In (fT^ the function F{t) describes the decay of polarization of a dielectric sample 
with time after sudden removal of steady polarizing electric field. The frequency- dependent 



dielectric constants of broad class of materials inc 
be interpreted in terms of the Williams- Watts 



uding polymeric systems and glasses may 



F^{t) = exp[-{t/Tr] 



31| polarization decay function: 

< (T < 1, 



(13) 



where exponent a and the time constant r depending on the material and fixed external 
conditions such as temperature T and pressure. 

The relation p2|l can be generalized for case of external field as we are interested in the 
time scale lesser than the time interval during which an external standing electrical wave 
may be considered stationary. After substitution e{uj) e(i^,g) and eg — > est{g) the Eq. 
()12|) we can transform and write in the form: 
e(u;,g) - e„ 



est(g) 



g[a,uj) 



g[a,uj) = -a 



X = ujT, s = t/r. (14) 



For example in Debye's classical theory of dielectric relaxation a 
the right side of (fT^ has the form (1 + iujT)~^. 

III. AVERAGE INTERACTION POTENTIAL BETWEEN STERIC SPIN 

CHAIN AND EXTERNAL FIELD 

As we have already shown in order to take into account the contribution of orientation 
effects into the polarization one has to calculate the total average value of the interaction 
potential of steric spin chain with the external field —{5U l^{Eq,VL, ipo)) (^L^^'\^y 

Taking into account the fact that the external field is low i.e. \E{lx — x)\ <^ \Eioc{lx~x)\ — 
\Eint{lx — we can apply the Taylor decomposition of dipole angular momentum: 

-x) -x) + -x), Sp{x) E{x;Eo,n,Lpo), (15) 

where |5p(x)| <^ |p°(x)|, as well as E^^^(x — l^) and — l^) are correspondingly the field 
strength and the dipole angular momentum of the molecule, located in the /^-th cell in the 
absence of the external field. In case when the coordinate x is outside of Ix-th cell, the field 
vanishes. Inside the cell they have constant values. From the discussion above and taking 
into account (|7j), the relation ()15p can be represented in the form: 

Lx Lx 

-6UL^{E\x,g)) = J2p''i^--^)E{x;Eo,n,^o) + J2Sp{lx-x)E{x;Eo,n,^o)- (16) 

Using the fact that without an external field the spin system has no polarization we conclude 
that first sum in p6|) vanishes. In other words, the interaction potential has the following 
form: 

-6Ul^ (E|x, g) = J2 ^P(^- - ^)^(^; ^0, ^, <^o). (17) 

Now we turn to the equation of motion for the steric spin-chain with relaxation in 3D spin 
lattice in external field. Recall that the interaction potential —6Ul^(Ei\x, g) (see ((Tj) and (fTTj) ) 
between ID disordered spin chain and external field does not take into account relaxation 
with environmental spin chains. It is possible only after solution of the dynamical problem. 
The resulting interaction potential in this case will be complex, where the imaginary part 
characterizes the relaxation processes in 3D lattice. 



1 Q, 
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the integral in 
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In more general case those motions can have quantum characters and can be represented 
by complex Langevin-Schrodinger type stochastic differential equation: 

X5UL^iE\x,g) = \E + <i/-\dt)^^, (18) 

where 

rriQ and fi are the molecule (spin) mass and the spin chain's effective mass correspond- 
ingly, t denote of natural parameter of evolution along the spin chain. In the equation 
fjlSp interaction potential f/i^(E|x, g) is a random complex function. Below we present its 
detailed description. 

Substituting 

^(t) =exp( /" E{t')dA, (19) 







into (|TH|l and using the relations we obtain the following nonlinear complex SDE Q]: 
Si + S2 + A(E-F) +A7(t) = 0, E{t) = 9{t)+i^{t), Et = dE/dt, (20) 

where 

Y,P{lx - x)E (x; Eo, n, (^o) = F + 7(t), V=(J2 ^Pi^- " ^)E{x; ^o, ^, <^o))^. (21) 

lx=0 lx=0 

In formulas (^0)) and (f7T|) we have denoted by V the mean value of the sum, and by f{t) its 
complex random part. Analyzing the contribution from different mechanisms of molecule 
polarization in glass medium, we conclude that under the influence of the external field 
Q with frequency Q ~ 10^ Hz the main part comes from elastic dipole (dipole thermal 
)olarization is not essential in this case due to the large relaxation time r ~ 10"^ ^ 10"^sec 
26^ . 3^). Let us note that the coefficient of elastic dipole polarization at low external fields 



is determined by 

^dipiL -x) = A"^ [p" sin[/?(/^ - x)]) = sin^[/5(/:j; - x)], 

Sp{l^ -x) = adipilx - x)E (x; Eq, fi, (/pq), (22) 

where j3{lx — x) is the angle between the external E (x; i?o, f^, V^o) and the internal E^^^{lx — x) 
fields, A(/^) = p°{QE-^^{Q = p^E^^^ is the dipole energy in the field Following 
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heuristic argumentation of Debby j26|, |3ll . l32l | one can write down the expression for the 



elastic dipole polarization, which takes into account the spin (polar molecule) relaxation 
process in the glass: 

i-ch{lx -x) = adipik - x)/{l - iUt) = ^0 _ ^Q^^ ' (23) 

where r is the spin relaxation time in the glass. It is very small in the aforementioned media 
r ~ 10^^^ ^ 10"^^ sec (see, for example (2(t||). 

So, the equation (|TH|) with generalized coefficient of elastic dipole polarizabilitiy will 
describe the motion of ID random steric spin chain in the external field with relaxation. 

After substitution of ()23|) in ()2ip and simple calculations we obtain: 

— , ^ -9. x\ 1 + N^p^'eI f sin(2fcLJ\ , , 

V = {j:r,,(L-.)E^.;EM)^ = -^^^^ (l - , (24) 

where Eq{(Po) = 2Eosirnpo. 

Let us now investigate the properties of the random function f{t). From the relations 
fETjl and (PHjl it is easy to find the random strength: 

lit) = -iV^ ■ ^(^) = + ^°<2fc*t)) p cos2^(/, - t) 



P{h-t)=P{k-x), kt = h~'^/2jlk, k = hl/^. (25) 

If the phase (3 is homogeneously distributed along the interval [0, it], then for the mean value 
we obtain: 

MJit) = (fit)) = 0, Mm = m) = 0. (26) 
or the autocorrelation function one can write the following expression (see, for example 



37|): 

2 / , . n— 2 \ 2 



- (7(07(0) ^ i (^) (Iflf) (Mm) . (27) 

where t = t — t' . Substituting in (P7j) and carrying out straightforward calculations, we 
may obtain: 
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So, we have shown that the dispersion equation (|T^. which describes the relation between the 
external field and the local polarization in disordered system, is reduced to the investigation 
of the Langevin type nonlinear complex stochastic SDE with the autocorrelation function 
fl26|l - PH|l of white noise. For further investigation, it is convenient to represent the complex 
equation as a system of two real equations: 

e + e'^ -^^ + X{E-Vr + frit)) =0, (29) 



^ + 2M + X{V, + Mt))=0, (30) 

where 6 = dtO, ^ = dt^, K = ReV, Vi = ImV, frit) = Refit), and /^(t) = ImJit). 
Now the problem is to find of evolution equation for the conditional probability: 

g(^,^,t|^o,^o,to) = {5{e{t)-e{t,))5{m-^{h))) , (3i) 

{6lo=e(to);i9o='f(to)} 

describing the probability of fact that, the trajectory (6*0 = 6'(to),'*9o = ''^('^o)) at the initial 
moment leaving of natural parameter to from the point (^oi'^^o); at the arbitrary moment of 
t it will be turned out to the suburb of point (6,^). Subject to SDE system (P^ . (jHUj) the 
Fokker-Plank equation can be easily found (see also 

f = + + - + A(E - + (2«tf + At/.)§ + m, (32) 

where Q = (5(E|6', t), Dr = X^ReD and Di = X^ImD. Note, that the solution of equation 
()32|1 must satisfy the initial condition: 



Qme,d-t) 



= 5{9-e,)5{§-do), (33) 

t=to 



where the initial phases 9q and are equal to zero. Later we will be more interested in 
the fixed limit of the solution ()32|1 . which is obviously received under values t ^ At = 0(1) 
which is equivalent to the condition t — > oo. In that limit the equation ()32|1 is simplified 
and accurate to the value 0{At/tN^) ^ 0(1) may take the following form: 

Dr^ + A^Q. + (^^ - + A(E - K)) ^ + {29^ + XV.) ^ + 49Q. = 0, (34) 

where Qs = Qs(E\9,-d) = limj^oo ^(6*, ^9, i) and t^v^ = A^^. 

The ()34|) is an elliptic type differential equation for which there are not real characteristics. 
This implies that the type of quasilinear equation depends on which solution is considered 
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and can be different for different solutions. The solution of the equation ()34|) must satisfy 
the conditions: 

0, \n\ = {e' + ^'f\ (35) 



OQs 



s dn 



s 

where the border condition is set on a curve S and n is correspondingly the normal of curve. 
Now one can write the value of the average interaction potential: 

1 1 Z"*'^^ / • • \ 

{6Ul,(E\x, g))^ = E + - ■ lim - / (9 + i^ + 9^ -^^ - i29^) dt 

A t^tN^ i Jo ^ ' 

= E + V ■ lim - \ (9 + 1^ + 9'^ -d"^ - i29A dt + f — ^ . (36) 
A t-^+^tJo V / \tNj 



Using the Birgoff ergotic hypothesis j30[ we can change integration by natural parameter t 
to integration by stationary distribution Qs(E\9, ??) and finally find the following expression 
for mean value of the potential energy increment: 

(5f/L,(E|x,g)),. = E+— y [9^ -^^ -t29^]QsiE\9,^)d9d^ + 0{N-'), (37) 
where R = Qs(Ei\9,'d)d9d'd denoted normalization constant. 

IV. STATISTICAL MECHANICS OF THE STERIC ID RANDOM SPIN CHAIN 

Before turning to the effective statical dielectric constant calculation the two important 
problems have to be solved. The first one is the calculation of the partition function Z(E; g) 
for the steric ID random spin chain with energy E, whereas the second one is the determi- 
nation of the thickness of layers on a scale of this chain (i.e. the fragmentation of the media 
on layers with different dielectric constant). 

First, let us consider the statistical properties of the ID steric spin chain. As it is shown, 
the type of a compound {a — Si02) can be used as a canonical rnodel for such network glasses 

(see Fig 1). It means that 



describing continuous random networks of atom and bonds l38 

Heisenbe., spin-gla. Ha^iltonian Q ca„ be ...table fo. desCbing of SO disordered spin 

system. Note that this allows to use the ergodic hypothesis (see ((Tj)) and make reduction of 
3D spin lattice dynamical problem to the problem of ID complex Langevin-Schrodinger type 
stochastic differential equation (see (fTH|) . ((201)) and the problem of construction statistical 
mechanics of steric disordered spin chain system. As it is easy to see, this type of reduction 
in a certain sense is conditional because both problems are interrelated. 
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For further investigations it is favorable to consider the spin-glass Hamiltonian of type: 

^0 

H(iV,; g) = - 5^ Mrr,)SiS, +p''Y. ^'^^^ = ^' (38) 

ij i ^ 

where rij = \i — j|(io(T) + rjij is a distance between ID spins Si and Sj (classical vectors 
of unit length), rjij is its random part, subject to Gaussian distribution with zero mean 
value and unit variance, Jjj is nearest-neighbor interaction constant, which depends on the 
distance between spins. It may be positive or negative, the external field Ei is defined by 
means of formulas (j21) • 

Now the problem is the construction of the energy distribution function of spin-chains 
Z(E; g). However, the problem is that the time scale, on which we make statistical study of 
the system, is very short (<10~^°sec = 0, Insec), while the characteristic thermal relaxation 
time in amorphous media j26| is of order Vt^^ ~ 10~^ -^ 10~^sec, where Vt^ is the frequency 
of thermal fluctuations. The last fact means that the temperature and related thermody- 
namical constructions become meaningless in our problem. Nevertheless, some structural 
similarity between gas and amorphous media is evident. At that a steric spin chain will 
correspond to an atom in gas. Since in the equilibrium state the average energy value per 
atom in the gas is f/cT, a corresponding value in this case will be the chain's energy in the 
equilibrium state (chain's energy without external field). However, the considered system 
has a specific peculiarity. The point is that the equilibrium state in gas is characterized by 
one temperature whereas the spin system can be in the equilibrium state in any negative 
energy. These energies coincide with the local minimums of non-perturbed Hamiltonian 



In other words, in this case 



and as it is well-known their numbers can be high 
the phase space can be decomposed uniquely into micro-canonical states associated with 
different thermodynamic equilibrium states j4l|. 

The Hamiltonian (pH|l in case of the external field absence can be rewritten in spherical 
coordinate system as follows: 

Ho(do(T),iVa.) = ^ Jij(rij)(|cos?/^jCos?/^j cos(</)i - </)j) + sin?/^jsin?/^j^, (39) 
{ij}=i, i^j 

For the determination of local minimums one has to solve the following algebraic equations: 
(0) = —— = ^ Jijirij) i - sin ilJi cos cos(</)i - (pj) + cos ipi sin iljj j = 0, 
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dH 



^r,, (® ) = = ^ (^cos ipi COS cos(0i - 0j) + siii ipi siii j =0, (40) 



where 0j = {ipi,(f)i) are angles of i-th spin (t/^j is the polar and 0j azimuthal angles), = 
(Oi, 92....0Ar^) correspondingly describe the angular part of spin configuration. Now suppose 
that non-perturbed Hamiltonian for fixed averaged distance between spins doiT) has n local 
minimums of function E((io(T)), each corresponds to Mj spin configurations {0°*}, where 
i = 0, l...M^. The number of all configurations, which correspond to local equilibrium states 



Mfuii = Yl!j=i Correspondingly, the weight of every equilibrium state may be defined 



by formulas: 



N 



P,(E,;do(T)) = MjMf^a, J] P.(E,; do(T)) = 1. (41) 



Thus we can propose to use the following statistical weight instead of the canonical 
distribution at multi-equilibrium state: 

H^„(H;g) = i'„(E.<<.(T))exp{-|(|^}. (42) 

where 'Ei^{dQ{T)) is the energy of spin chain in the absence of an external field. Recall, that 
the multicanonical ensemble was introduced in the |3| as an approach to simulate a strong 
first-order phase transitions. It has been very successful in sampling a large barrier. 

Now, taking into account (j39p - (j42|) . one can present the expression of the energy distri- 
bution function for spin system (this expression can be explained as a local partition 
function in above mentioned multicanonical thermodynamics): 

Z„,E„;,„).P(E.4(r,) /§...;5^e.p{-|i^}. (43) 

where q describes the set of random distances rjij and random angles (3i, correspondingly 
dVLi is an element of the solid angle f2j containing the unit vector S^. 

Energy distribution function of spin chain can be essentially simplified after the suggestion 
that in Hamiltonian (jHH|) only nearest-neighboring spins interact, i.e. Jij = if |j — > 2. 
In this case the multidimensional integral can be taken exactly as follows. The integration 
starts from the end of the chain. When integrating over dVL^ we take as a polar axis the 
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direction of the vector {Jiu_i^Si + p^Ei). Then it is easy to obtain the following expression: 



Z^(E^;g,q) = P^TT - / exp{i^i;^cos^}si 
i=i 



sin %l) dip 



Nx 



i=l 



E. L 



Jl^., + 2p"E,J(,_i),cos A + 



(44) 



where (3i is the random angle between the vectors Si^i and Ei. Supposing that the distri- 
bution of spin Si around field Ei direction is isotropic, one can perform an integration by 
angle After simple integration with the help of formula j^]: 

Ei(ax) = I — dx = In |x| + > — — — , a 7^ 0, 



k=l 



we can find: 



Nx 



Z(E;g,r/)=P(E;do(T))n 



i=l 



1 r sinh Ki 



2 Jo 



sin Pi dPi 



Nx 



P(E;do(T))n 



i=l 



Aai 



Ei{hi + Oj) — Ei{—bi — a. 
2p°£'jJ(j_i)j, 



Ei{bi — tti) — Ei{—bi + Oj 



(45) 



For simplification of expression and below will be admit index z/. In the fl45|l symbol rj 
describes the set of random distances, Ei{x) is the exponential integral function. It is simple 
that in spite of the fact that every member in parentheses near the points E^ = Jo/p^ or 
E~ = —Jo/p^ has a singularity nevertheless their sum is an analytical function. 



In the paper |4^ at the first time it was suggested that one may describe spin glass by 
a Hamiltonian of the type (jHHj) . where spins are put onto the sites of a regular lattice, and 



disorder is introduced by a suitable distribution W(J(i-i)i) of exchange bonds. A standard 
choice is the Gaussian Edwards-Anderson model kjl (see also 



I^(J(i-i)^ 



27r(AJ( 



: exp 



{J{i-l)i ~ Jq) 

2(AJ(,_i),)2 



Jo - {J{i~l)t) (AJ(,_i)i) - (4-l)*>at. " {'^i^-^)i)av- 



(46) 



Recall that for this model Jo and AJ(j_i)j are independent on distance and scaled with spin 



number as 



(J(,_i).>^„ = Jo cx N-\ AJ(,_i), oc N-'l\ (47) 
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to ensure a sensible thermodynamic limit. Eq.s ()46p and ()47|) (•••)^^ describe the averaging 
procedure. Now we can conduct averaging of function ()45|) by distribution pUj) and find the 
local partition function near some equilibrium energy of disordered spin chain: 

Z(E;g) = P(E;rfo(T))(Z(E, J;g))^, 

A': 



G(E,J,;g) 



2a,; 



(Z(E,J;g))^ = f] f rfJ(,_i),H^(j(,_i)i)G(E,J,;g), 
i=i ^ 

Ei{bi + ai) - Ei{-bi - flj) - Ei{bi - a,/) - Ei{-bi + a, 



(48) 



Taking into account the fact that on the scale of half-wavelength A'j. ^ 1 as well as that 
G'(E, Jj; g) is an analytical function for computation of integrals in partial partition function 
the Laplace asymptotic method may be used: 



Z(E, g) ^ P(E; doiT)) n EiilP, + a°) - Ei{-b\ 



i=l 



^° a?) 



Ez(60-a°)-Ez(-6^ + a^) 



I 



i2/£;.Jo, 6° = i 



(49) 



It is simple to show that in the limit of -Ej — > the correspondence Z(E, g) -P(E; (io(^)) 
takes place. 

Note, that averaging procedure in the framework of multicanonical ensemble is carried 
by the following formula: 

, E:Li>^(E.|x,g)Z(E.;g) /ly(E|x,g)Z(E;g)c^E 



Z(E; g)dE 



E:=i Z(E.; g) 

where Y{x, g) is an averaged physical parameter. 

For ID random steric spin chain system in analogy to usual thermodynamics Helmholtz 
type free energy may be specified near the local equilibrium state with energy Ej,. In this 
case the free energy amount with the expectation of one spin and the per unit of equilibrium 
state energy 'Eiy is defined: 

1 



F(E;g) 



lnZ(E;g) 



(51) 



All macroscopic thermodynamic properties of 1-D random steric spin chain can be obtained 
by free energy derivatives. The simple investigation of expression for free energy shows that 
in the steric ID random spin chain model with the nearest-neighbor interaction depending 
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on amplitude of external field £"0 there are not phase transition phenomena. However in this 
case under the even low fields Ei ~ l/N^ the free energy of spin system essentially changes 
(see Fig 6). 




FIG. 6: The Helmholtz free energy for any average energy of spin chain (see figures (a),(b) 
and (c) ) and amplitude of external field £"0. On the figures axis k shows the parameter which is 
proportional to external fields amplitude Eq. 

Particularly, the order parameter for disordered spin-glass medium is described by ^ ■ 
and without the external field on the scale of the period of standing wave it is equal to zero, 
where is a projection of spin on the direction of external field propagation. It becomes 
non-zero at the weak external field ~ \E^\ due to a symmetry breaking. In this case 
nanoparticles (spin chains) with the super spins (macropolarizations) are generated on the 
microscale space and nanoscale time. 

V. DIELECTRIC PERMITTIVITY OF THE NEIGHBORING LAYERS 

Now we can examine the question of dividing the dielectric medium on the regions with 
extremely different polarizations. Recall that the polarized medium without an external 
field has a zero macroscopic polarization on the space scale 10~^ ^ 10~^cm. After turning 
on the external field in the medium (on the scale of wavelength of standing electrical field) 
the amplification of the orientational effects and, as a result, the initiation of macroscopic 
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polarization is possible. Assuming that the collective orientation effects, generated in those 
regions of media, where voltage of external field E is of order to the critical value E^, we 
can divide the scale of wavelength of external field on the forth region. In the first and third 
regions the macroscopic polarization of medium is zero whereas in the second and forth 
regions it is different from zero (see Fig. 7). Note that the height of these layers may be 
computed with the help of numerical experiments. 

E, 






FIG. 7: The regions with different polarizations and correspondingly different dielectric constants 
on the one wavelength As scale. 

Now, with calculation of effective dielectric constant at the half- wavelength scale of ex- 
ternal electromagnetic field at the time intervals At ^ 2ti/VL. Using (fTT|) . we can obtain: 

Remember, that here we have taken into account the contribution from the elastic dipole po- 
larization only. It is easy to see that the nominator of the expression (|52p has no singularities 
and is bounded by the imaginary part of the potential increment ()37p. The important char- 
acteristic is the dielectric constant's difference in neighboring layers. Taking into account 
the relations (HH) and ^ one can get: 

^^..(g) = 6.,(g) - 6s,(0) = (1 - 6s,(0)) + Y^^- (53) 

Dielectric permittivity of layer in the external field can be calculated simply: 

est(g, cu) = Coo + (erf(g) - Coo) Q{(y, i^)- (54) 
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From equation ()54|) taking into account Debye relaxation model in the x-rays region we find: 

6.,(g,^) ^ 1 - (6,,(g) - + (^)'}' (55) 

where ujp is the plasma frequency of media. 

Now we can write the expression for difference of frequency-dependent dielectric constants 
(dielectric permittivity function) in neighboring layers: 

5e(g, u) = estis, uj) - est{uj) = 5e,t(g) ^'(cr, uj), (56) 

which in the x-rays region has a form: 

56(g,^)^-5e.,(g)|.^+(^)'}. (57) 

So, we constructed the dielectric permittivity function for spin-glass medium under 

the external standing electromagnetic wave (0), which is correct on the micrometer space 
scale and nanosecond time scale. 

Note, that the main peculiarity of expression (jSlj) concludes in the fact that est(g) is a 
complex function. Moreover depending on parameters of weak external field est{g) can have 
big imaginary value. This fact in contrast to usual case builds the media with radically 
new properties. Recall when the external field is absent the dielectric permittivity function 
depends on e^, which is order of unit real positive constant. 

VI. CONCLUSION 

In this article a new microscopic approach is developed for the study of the properties of 
stationary dielectric constant and dielectric permittivity function in the dielectric medium 
under the external standing electromagnetic field. The approach consists of the following 
two general steps: 

1. Generalization of Clausius-Mossotti equation for dielectric constant in the external 
standing electromagnetic wave; 

2. generalization of the equation for dielectric permittivity function taking into account 
the previous results. 
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Mathematically the problem is solved as follows. The dielectric medium in the external 
electromagnetic field is modelled as a 3-D spin-glass system under the external field. Note 
that all general changes of properties of media take place in the scale of wavelength of 
external field. For that the layer of medium which consist from disordered ID spin chains 
with length of order of wavelength in detail has been investigated. Taking into account the 
fact that the distribution of spin chains in the infinite (x, y) plan is isotropic we can use the 
Birgoff ergodic hypothesis (see ©-©) and conditionally reduce the initial 3-D spin-glass 
problem to two ID problems. It means that we can investigate two ID problems separately 
only in this case the parameters of the first ID problem must be taken into account during 
the solution of the second ID problem. 

In the work all formal definitions for accounting the investment of orientational effects in 
the stationary and frequency-depended dielectric constants computation are adduced. 

The first ID problem is related to one-dimensional disordered quantum A^-particle system 
with relaxation. The investigation of the motion in this system takes place in the framework 
of complex Langevin-Schrodinger type SDE (fTH|l . which can be transformed to 2D system 
of nonlinear Langevin type SDE ()29|) - (p^ . For probability distribution of interaction po- 
tential in ID spin chain with certain energies E and external standing electromagnetic wave 
Fokker-Plank equation ()32j) is obtained, which is found after using white noise model for 
stochastic forces and system of SDE (P^ - (|Hn|) . In the limit of long-range distance t oo, 
the probability distribution Q(E\6,{}]t) tends to the limit of its stationary limit Qs(Ei\6,i)) 
which satisfies the elliptic type differential equation ()34j) . Note, that this type of differential 
equation (jH^ depending on parameters can have extremely various solutions between which 
there aren't smooth transitions. In another words this type of behavior in the solutions we 
can interpret as a critical phenomenon in the 3-D spin-glass system. 

The second ID problem includes the computation of ID random spin chain energy dis- 
tribution in the external field. In other words, it means the development of multicanonical 
statistical mechanics for disordered ID spin chain system in the external field. In order to 
do this the classical spin-glass Heisenberg type Hamiltonian is investigated (jSHI)- In the first 
step, the nonperturbed Hamiltonian (jH^ is used and the system of algebraic equations (jlUj) 
is found for computation of all possible stable spin configurations on the spin chain scale. 
After this, the statistic weight P^CE,^] dolT)) for certain energies E^, of stable non-perturbed 
spin chain is found simply by formula With the help of formulas (jlT|) - (PH|l the partition 



23 



Medium 



Vacuum 





Vacuum 



FIG. 8: The regions of different polarizations. 



function Z(E; g), which gives an energy distribution in the spin chain system after inserting 
of external field, is defined ()48|). The Helmholtz free energy ()5H) is constructed and it is 
shown that the ID random spin chain system subjected to weak external field Ei and the 
order ccl/N^ can be severely changed (see Fig. 6). 

The mean value of complex interaction potential (5?7l^(x, g))(|) between spin chain and 
external field is found from formula (0) taking into account expressions (jlT|) and PH|) . An 
averaging potential is used for computation of polarizability coefficient related with the 
orientation effects and correspondingly for generalization of Clausius-Mossotti equation 
flll|) . Remind that Eq. (fTTj) has a content only on the micro-scale space and the nano-scale 
time. Note, that the analysis of Eq. together with (f^Hj) show that in the spin-glass 

medium with the help of weak external field it is possible to form regions (layers) with 
different stationary dielectric constants est(g) (see Fig 8). As it is easy to understand the 
value of stationary dielectric constant strongly depends on initial electrostatic dielectric 
constant of media. In particular, it is shown that in the spin-glass mediums with static 
dielectric constants < 4 under the external field the value of the stationary dielectric 
constant can be changed maximum 3-4 times. In the case > 4 analysis shows that the 
value of -Re(A(g)) 1 and correspondingly stationary dielectric constant can be complex 
function esj(g) depending on the parameters of external field g. Moreover for a set of some 
parameters of weak external field it is safe to say that the est(g) can have a large imaginary 
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part (see 

Obviously, near to this value critical effecters are hold. In this case taking into account 
the relaxation processes, which are going in the lattice, the system becomes ordered and is 
characterized by macroscopic classical polarization. 

The second main step of investigation is the frequency-dependant dielectric constant or 
dielectric permittivity. Generalized Clausius-Mossotti Eq. (fTT|) may be used for derivation 
of generalized equation for dielectric permittivity ()14|) for the frequencies range u ^ Q. 
It is important to say that the analysis of stationary dielectric constant value est{s) shows 
that dielectric permittivity function e{w, g) (see (jSH) and (jHKjl ) can have positive, as well 
as negative dispersions. The last circumstance is very important especially for the x-range 
frequency ui because on the scale of wavelength of external field (see Fig. 7) the dielectric 
permittivity in the region of positive dispersion "I" and "III" Reei{uj, g) < 1 whereas in the 
neighbor re gion of negative dispersion "II" and "IV" correspondingly Ree2{uj,g) > 1 (see 
for example |4^ 0, 1^ ) • Recall that this occurs when initial electrostatic dielectric constant 
es > 4. Evidently that the difference between dielectric permittivity of neighboring layers 
(|K7jl in the x-ray region as compared with usual cases when the external field is absent can 
differ drastically. 

In this case collective orientation effects of spins make main investment in the permittivity 
function (macroscopic polarization) which radically changes behavior of media on the space- 
time scale of external standing electromagnetic wave. Particularly only a presence of large 
imaginary part in the coefficient of Sestis) brings to an essential growth of the real part of 
difference {^i^Up/uj) + {up/ujYY 

Thus, obviously it is possible to make and control the periodic structure from dielectric 
permittivities in the spin-glass type dielectric medium even with the help of weak external 
standing Q electromagnetic field. In other words a new type of controlled radiator for 
generation of coherent x-radiation is suggested. 

On the other hand the periodic dielectric permittivity function may be illustrated as a 1-D 
lattice of quantum dots (quantum strings) which can be interesting in the view of quantum 
computer making. 

Finally it is important to note that the first time a new mathematical method for re- 
ducing the problem of dynamical 3-D disordered spin system in the external field to the 
two conditionally separated ID problems is developed. The last circumstance can be used 
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for elaboration of new extremely effective parallel algorithms which is very important for 
systematic investigations of mentioned problem in the framework of numerical simulations 
method. 
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